Energy carriers in the Fermi-Pasta-Ulam (3 lattice: Solitons or Phonons? 
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We investigate anomalous energy transport processes in the Fermi-Pasta-Ulam f3 lattice. They are 
determined by the maximum sound velocity of the relevant weakly damped energy carriers. That 
velocity can be numerically resolved by measuring the propagating fronts of the correlation functions 
of energy/momentum fluctuations at different times. The numerical results are compared with the 
predictions for solitons and effective (renormalized) phonons, respectively. Excellent agreement has 
been found for the prediction of effective long wavelength phonons, giving strong evidence that the 
energy carriers should be effective phonons rather than solitons. 
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Energy transport in low dimensional systems has at- 
tracted enduring interest [J-Q- O ne striking finding is 
the phenomenon of anomalous transport [5|-|7[, which 
has recently been experimentally verified for carbon nan- 
otubes [|[ . Theoretical efforts [lllliol usually follow the 
pioneering work of Peierls and focus on the low tempera- 
ture region, where weakly interacting phonons are consid- 
ered to be the responsible energy carriers. A microscopic 
transport theory beyond the low temperature regime is 
still lacking, which leaves the explanation for most ex- 
isting numerical and experimental results far from sat- 
isfactory. One central question concerns the type of en- 
ergy carriers whose properties determine the underlying 
transport behavior at higher temperatures. As tempera- 
ture or nonlinearity is increased, collective motions other 
than phonons could also be excited. It is thus desirable 
to identify the specific energy carriers for these low di- 
mensional systems. 

The Fermi-Pasta-Ulam j3 (FPU-/?) lattice is a classic 
example showing the effect of anomalous transport, and 
therefore a perfect testbed for comparison between the- 
oretical predictions and numerical experiments. Anoma- 
lous transport manifests through an increase of the heat 
conductivity with the system size. This in turn im- 
plies that the responsible energy carriers are anomalously 
weakly damped and propagate "ballistically" over very 
long distances. Due to nonlinearity, solitons fi"i"l.[l2l . dis- 



crete breathers [I3| and interacting phonons [1J-|18| are 
candidates for these carriers. In particular, supersonic 
solitons have been considered as major energy carriers 
which are responsible for the anomalous transport be- 
havior 0, [20[ . The pivotal evidence supporting the idea 
of soliton transport is the numerical observation of ul- 
trasonic energy transfer. The sound velocity c s of en- 
ergy transfer was measured by following the spreading 
of an initial energy pulse, using both non-equilibrium 
[l9l . l2l| and equilibrium [20[ methods. In particular, the 
temperature-dependent sound velocity c s is compared 
with a prediction derived from soliton theory [1 2j and 



good agreement has been found in Rcf. [T9[. However, 
strong finite size effects of the soliton velocities were not 
clarified in Ref. [L2\ . The same data for c s in Ref. [lj| 
are also in good agreement with the predicted velocity for 
effective phonons fl7| . The uncertainty of the computed 
data in Ref. [l|| is too large to distinguish between the 
two predictions. To identify the true energy carriers, a 
more accurate numerical determination of c s is needed. 

In the present paper, we apply the equilibrium ap- 
proach recently developed by Zhao [2Cj to study the 
energy transport properties in the FPU-/? lattice. As 
demonstrated in Ref. [20| , the sound velocities c s can be 
measured with very high precision. We will show that 
the numerical results are in very good agreement with 
the prediction for effective phonons and the agreement is 
not limited to the FPU-/3 lattice. The soliton predictions 
show clear deviations and can be ruled out. 

We consider the dimcnsionlcss Hamiltonian for the 
FPU-/? lattice 
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where pi denotes the momentum and qi denotes the dis- 
placement from equilibrium position for the i-th atom 
with i = 0,±1,±2, ...,±N/2. The local energy den- 
sity is defined as Hi = pf/2 + — qi-i) 2 /2 + (<& — 
<7i_i) 4 /4. Fixed boundary conditions arc applied to the 
two end atoms, which are additionally coupled to stochas- 
tic Langevin heat baths with specified temperature T. 
The normalized correlation functions of energy and mo- 
mentum fluctuations are defined as uOf 



C E (i,t) = 
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where Ai?*(t) = Hi(t) — (Hi), and (•) denotes ensemble 
averages. Note that CE/p{i,t = 0) = 5^o. Therefore, the 
correlation functions C^/piht) describe the spatiotem- 
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FIG. 1: (color online) Spatial distributions of the correlation 
functions of Cs(«,t) (upper panel) and Cp(i,i) (lower panel) 
at t = 60 for T = 0.02, 0.5 and 5 for the FPU-/3 lattice with 
TV = 500. 



poral spreading of the initial energy/momentum fluctua- 
tions [20| . 

In Fig. [TJ we depict the spatial distributions of the 
correlation functions CE/p(i,t) at t = 60 for different 
temperatures. The values T = 0.02, 0.5 and 5 correspond 
to low, intermediate and high temperature regimes, ac- 
cording to the scaling behavior of the aspect ratio e = 
((<& — 9i-i) 4 ) / (fe — Qi-l) 2 )- The aspect ratio scales as 
e oc T in the low temperature regime and as e a y 1 / 2 
in the high temperature regime. As can be seen from 
Fig. [TJ both distributions possess symmetric propagat- 
ing fronts with identical propagation velocities. These 
propagating fronts are induced by the fastest travelling 
energy carriers 0, [2(J . With increasing temperature the 
propagation velocity increases, which is caused entirely 
by the presence of nonlinear terms in the equations of 
motion. The fluctuations of Cp(i,t) are much smaller 
than those of CE(i,t). Therefore, we will determine the 
sound velocity c s by measuring the peak positions of the 
propagating fronts for Cp(i,t), as in Ref. [2(j. We have 
also tested that the results do not depend on the type of 
boundary conditions (fixed or periodic). 

Let us briefly introduce the predictions of sound ve- 
locity for solitons and effective phonons. According to 
[IHIl2j|. the soliton profile with z ~ — is given by 



Q 



,(z) = v/2( C 2 - l)sech (2 ZV /( C 2 - l)/ci) . (3) 



It f ollows t hat the energy of a soliton is proportional to 
c s \/ c s — 1 1 where we assume c s > without loss of gener- 
ality. Using a Boltzmann distribution for energies of ex- 
citations, we conclude that the temperature T introduces 
an energy scale such that larger soliton energies are expo- 
nentially suppressed. At the same time smaller energies 
imply smaller sound velocities. Therefore a Boltzmann 
distributed gas of solitons will typically show maximum 
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FIG. 2: (color online) Sound velocity c s as the function of 
temperature T for the FPU-/3 lattice with N = 1000. The 
circles are the numerical results. Error bars have been plotted, 
but they are much smaller than the symbol size. The dotted 
line is the analytic prediction for effective phonons from Eq. 
([5]) and the dashed line is the analytic prediction for solitons 
from Eq. ([4]) with ri = 2.215. The numerical result obtained 
for T = 0.5 in Ref. [20] is plotted with a cross symbol. 



sound velocities which correspond to solitons with an en- 
ergy of the order of the temperature T. Therefore it 
follows that [H 



(4) 



The constant a is a free fit parameter. Note that c s — > 1 
as T -> and c s w ? ? 1 /4 r i/4 as y » 1 in the high 
temperature regime. 

For the effective phonons, the sound velocity is de- 
fined as the maximum group velocity of the renormal- 
ized phonons. Here, rcnormalization implies a mean field 
treatment of nonlinear terms in the equations of motion. 
As a result eigenfrequencies of phonons are renormalized, 
and will increase with increasing temperature. There- 
fore, renormalized phonons will also yield sound velocities 
which increase with increasing temperature, becoming 
supersonic as compared to the case of T — > 0. In partic- 
ular, c s — du)k/dk\k=o where Cju = 2y / asinfc/2 with < 
k< 2tt and a = 1 + (£<(?< - %_i) 4 ) / (£ 4 fe - q t -i) 2 ) 
[l7|. It follows [17| 
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For T — > we find c 
ture region c, 



1, and in the high tempera- 



yj r °° x^e-^/^dx/ r °° x^e-^/^dx w 

1.22T 1 / 4 . Both predictions (solitons and phonons) yield 
three similar results for the sound velocity: (i) c s > 1; 
(ii) c s (T — > 0) — > 1; (hi) in the high temperature regime, 
the sound velocities exhibit the same scaling with tem- 
perature as c s oc T 1 / 4 . 

In Fig. [21 we plot the numerically determined sound 
velocity c s as the function of temperature T for the FPU- 
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FIG. 3: (color online) Spatiotemporal evolution of energy 
densities Hi(t)/T (upper panels) and relative displacements 
qi(t) — qi-i(t) (lower panels) at thermal equilibrium. The left, 
middle and right columns correspond to the harmonic lattice 
at T = 1, and the FPU-/? lattice at T = 1 and T = 20, re- 
spectively. The lattice size N = 100 and periodic boundary 
conditions are applied. 

/? lattice. The computational errors are extremely small 
and c s is measured very accurately. The numerical results 
are compared with the predictions for solitons from Eq. 
((31), and for effective phonons from Eq. (JSJ). Excellent 
agreement has been observed for the effective phonon re- 
sult in the entire temperature region being explored. The 
soliton curve corresponds to r\ = 2.215, which reproduces 
the correct high temperature result. However, the devi- 
ation from the prediction for solitons of Eq. (J4j) is quite 
distinct in the intermediate temperature regime where 
0.05 < T < 2. We note that variations and optimizations 
of rj do not improve this discrepancy. The cross symbol in 
Fig. [2] represents the sound velocity c s = 1.31 measured 
at T = 0.5 in Ref. [2(|. Although it was originally viewed 
as an evidence for soliton transport, the numerical result 
is actually coinciding both with our numerical data and 
with the effective phonon results, in contrast to the soli- 
ton theory. This finding provides strong evidence that 
the effective phonons, rather than the solitons, should be 
the energy carriers responsible for anomalous transport 
in the FPU-/? lattice. 

To visualize the energy transport processes, in Fig. [3] 
we plot the spatiotemporal evolutions of local energy den- 
sities Hi(t)/T and the relative displacement qi(t)— <7i-i(i) 
for the harmonic lattice and the FPU-/? lattice, respec- 
tively. The systems are thermalized at a given temper- 
ature and then the heat bath is removed. The evolu- 
tion functions are recorded at thermal equilibrium and 
results for a suitable time window are displayed. The spa- 
tiotemporal patterns for both the FPU-/? lattice and the 
harmonic lattice are qualitatively similar. Since solitons 
are definitely excluded for the harmonic lattice, there 
appears to be no signature for soliton transport in the 
FPU-/? lattice as well. The qualitative behavior of the 



FIG. 4: (color online) Sound velocity c s as the function of 
temperature T for the H n models with n = 3, 4 and 5. The 
symbols correspond to the numerical results at N = 1000, 
whereas the lines are the predictions for effective phonons 
from Eq. (0. The errors are much smaller than the symbols. 

spatiotemporal evolutions does not change for larger time 
windows. 

In the high temperature limit, the FPU-/3 model can 
be reduced to an H n model with n — 4 assuming the 
following Hamiltonian 

N 

H n = Y i [p 2 i /2 + \q i -q i ^\ n /n\ . (6) 

i=l 

To demonstrate the power and consistency of the effec- 
tive phonon formulation, we consider three different cases 
with n = 3, 4 and 5. Following the same procedure as for 
FPU-/? lattice, the sound velocities of effective phonons 
can be expressed with a compact formula 

c s = [T((n + l)/n)/r(3/n)] 1/2 (nT) 1 ^/" (7) 

These predictions are plotted in Fig. U and compared 
with numerical results. We again find quantitative agree- 
ment for all three models. 

Let us discuss some details in the spatiotemporal de- 
pendence of the correlation functions. In Fig. [T] the 
distribution functions of C^/pii, t) show many peaks be- 
tween the two propagating fronts at T = 0.02. These 
peaks are typical for coherent phonon propagation in har- 
monic lattices. For T = 0.5 and 5, there are no visible 
additional peaks for both distributions and a big hump 
emerges in the interior region for CE(i,t). The disap- 
pearance of the intermediate peaks can be therefore at- 
tributed to relaxation processes. Indeed, the distribu- 
tion functions CE/p{i,t) for T = 0.5 at various spread- 
ing times t = 10,30,50 and 70 in Fig. [5] show that the 
short time behavior at T = 0.5 is similar to that for 
T = 0.02. The phonon modes relax faster at higher 
temperatures, while on the other hand, the long wave- 
length phonon modes possess very long correlation times 
even at very high temperatures. According to Ref. [20j |. 
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FIG. 5: (color online) Spatial distributions of the correlation 
functions of Cs(«,t) (upper panel) and Cp(i,i) (lower panel) 
at T = 0.5 for various values of time t = 10, 30, 50 and 70, for 
the FPU-/3 lattice with size N = 500. 

the correlation function of energy fluctuations C^(i,t) 
is nothing but the energy density probability distribu- 
tion function (PDF). Therefore we can study the en- 
ergy diffusion process by measuring the mean square dis- 
placement (MSD) as (r 2 (t)) = Y, t i 2 C E {i,t). For PDFs 
shown in Fig. [5j the MSD is dominated by the area 
around the propagating fronts and can be approximated 
as (r 2 (i)) oc t 2 ~ v where the exponent v characterizes 
the dcminishing of the peak area with time as a conse- 
quence of slow but unavoidable dcphasing of even long 
wave length pohonons. According to [2(| the exponent 



< v < 1, and therefore the diffusion process is su- 
perdiffusive as (V 2 (i)) oc t a with a = 2 — v. It is interest- 
ing to analyze the connection between superdiffusion and 
anomalous heat conduction for the FPU-/3 lattice. Recall 
that the heat flux of phonons is </& = v^E^ where k 
denotes the wave number. The correlation function of 
the total heat flux (J = Jfc) can be approximately 
obtained as (J(i)J(O)) « (c 2 s E (t)E (Q)) cx t~ v since the 
only long time correlation is due to the energy carried 
by the long wave-length phonon modes. Applying the 
Green-Kubo formula for heat conductivity [2| , we obtain 

K cx f* /c " (J(t)J(O)} dt cx N 13 where 0=1- v. Without 
knowing the exact value of (3 and a, we obtain that the 
superdiffusion and anomalous heat conduction are con- 
nected via the exponent relation (3 — a — 1 0, [22[ • 

In conclusion, we have investigated the energy trans- 
port processes for the FPU-/? lattice using an equilib- 
rium approach. To identify the energy carriers, we ac- 
curately measure the sound velocity of the energy carri- 
ers by following the correlated spreading of the initial 
energy /momentum fluctuations. The sound velocities 
are found to be in excellent agreement with theoretical 
predictions for effective phonons. This predicability has 
been further confirmed for a series of H n models. On the 
other hand, no signature of soliton transport has been 
detected by visualizing the spatiotemporal evolutions of 
local energy densities and relative displacements. There- 
fore our numerical results clearly reveal that the energy 
carriers are long wavelength phonons for the FPU-/3 lat- 
tice. 
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